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ABSTRACT 

This  paper  presents  and  demonstrates  the  enhancement  of  noise  tolerance  exhibited  by  the  improved  Max-Min  scanning 
method  (IMMS)  for  phase  determination.  Numerical  simulations  and  various  tests  indicate  that  the  IMMS  has  a  good  tolerance 
to  certain  ranges  of  noise.  One  of  the  IMMS  applications,  an  online  method  to  calibrate  the  phase  shifter,  is  presented  as  well. 
This  improved  algorithm  provides  an  easy  and  convenient  way  to  inspect  the  linearity  and  the  non-uniformity  of  the  phase 
shifter.  It  is  suitable  for  use  in  assessing  the  reliability  of  the  phase  shifting  setup  in  real  time. 

INTRODUCTION 

The  phase  shifting  technique  is  a  key  step  in  optical  measurement.  An  improved  Max-Min  scanning  method  has  been 
described  by  Ding  and  Cloud  [1 , 2].  The  principle  of  it  is  briefly  reviewed  in  the  next  section.  There  is  always  a  serious 
concern  about  the  environmental  tolerance  of  all  kinds  of  phase  determination  techniques,  because  the  noise  tolerance 
establishes  the  limitations  on  the  range  of  application  of  the  phase  shifting  techniques.  The  majority  of  current  techniques 
require  a  vibration-isolated  table.  This  paper  shows  that  the  application  of  low-pass  filtering  and  curve  fitting  techniques  to  the 
recorded  signals  in  the  time  domain  improves  the  noise  tolerance  of  the  developed  IMMS  method.  Numerical  simulations  and 
tests  are  presented. 

One  useful  application  of  the  IMMS  approach  is  presented  in  this  paper  also:  calibration  of  the  phase  shifter.  Generally,  a 
phase  shifter,  which  is  usually  a  mirror  mounted  on  a  piezoelectric  transducer  (PZT),  is  used  to  vary  one  beam  path  length 
artificially  so  as  to  vary  the  phase  difference  between  the  interfering  beams  by  some  supposedly  known  amount,  and  then 
some  intensity  image  signals  are  captured  to  calculate  the  expected  phase  information  [2, 3, 4, 5, 6, 7].  Although  some  of  these 
techniques  were  claimed  to  be  insensitive  to  phase  shifting  errors,  an  efficient  real-time  calibration  of  the  phase  shifter  is  still 
desired  to  understand  the  effectiveness  of  the  phase  shifter  and  to  make  phase  shifting  techniques  more  accurate.  The 
working  status  of  the  phase  shifter,  such  as  the  demanded  driving  voltage,  non-linearity,  stability  of  the  PZT,  and  the  tilt  of  the 
mirror,  will  more  or  less  affect  the  accuracy  of  the  measurement  results.  In  the  past  decades,  several  calibration  methods 
have  been  developed.  They  have  been  categorized  into  two  classes:  fringe  tracking  and  the  Carre  method  as  described  by 
Hedser  Van  Brug  [8,9].  The  fringe  locking  method  is  mainly  suited  with  an  accurate  2k  phase  step.  The  Carre  method  uses 
several  phase-stepped  images  to  calculate  the  phase  shifting  angle  by  assuming  that  the  phase  shifter  is  linear  and  uniform, 
which  is  likely  not  true.  Neither  of  these  two  methods  is  able  to  provide  the  linearity,  the  tilt,  and  the  non-uniformity  of  the 
phase  shifter. 

In  this  paper,  a  new  real-time  calibration  method  is  proposed  to  calibrate  the  phase  shifting  apparatus,  as  well  as  to  inspect  the 
non-linearity,  tilt  and  non-uniformity  of  the  phase  shifter's  movement. 

IMMS  METHOD 

The  principle  of  the  IMMS  method  is  briefly  given  as  follows  [1].  The  intensity  equation  of  an  interferogram  is  written  as 


I  =  Ix  +  I2  +2y[TJlCos  (</>) 


(1) 
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where  /  is  the  intensity  at  one  pixel,  /i  and  k  are  the  intensities  of  the  two  interfering  beams  respectively,  and 
difference  between  these  two  interfering  beams  at  that  pixel  [3]. 

From  equation  (1),  the  phase  angle  can  be  derived  as: 


'  is  the  phase 
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In  this  equation,  /  represents  the  intensity  at  one  pixel.  lMax  and  lMin  are  respectively  the  maximum  and  minimum  values  of  the 
pixel  intensity  at  the  same  pixel  position  as  /. 

By  changing  the  driving  voltage  to  the  PZT  phase  shifter,  the  phase  difference  ^  between  the  two  interfering  beams  can  be 
artificially  varied,  and  at  the  same  time,  the  digital  pictures  of  the  light  intensity  of  the  interfering  speckle  image  for  each  mirror 
location  are  recorded.  After  recording,  plot  the  light  intensity  values  versus  driving  voltage  at  each  pixel  position.  Low-pass 
filtering  and  least-square  curve  fitting  methods  are  used  to  remove  noise  and  fit  the  intensity  curve.  The  smooth  intensity 
variation  curve  is  obtained.  On  this  wave  curve,  the  gradient  is  computed  at  each  point,  and  then  the  peak  values  and  valley 
values  are  determined,  these  being  the  maximum  and  minimum  intensity  values  respectively.  The  phase  valued  at  each  point 
can  be  determined  through  equation  (2). 

On  the  fitting  curve,  the  derivative  value  at  the  point  where  driving  voltage  equals  zero  is  calculated  to  determine  the  quadrant 
of  <!>.  If  the  derivative  is  positive,  then  reverse  the  sign  of  phase  ^  obtained  from  equation  (2)  to  put  the  phase  angle  in 
quadrant  3  or  quadrant  4.  Otherwise,  put  the  phase  ^in  quadrant  1  or  quadrant  2.  The  above  processing  is  performed  at  each 
pixel  over  the  entire  image.  The  phase  map  is  determined. 


THE  NOISE  TOLERANCE  PROPERTY  OF  THE  IMMS  METHOD 

Since  some  numerical  filters,  for  instance  a  low  pass  filter,  can  be  employed  to  process  the  recorded  intensity  signals  in  the 
time  domain,  a  good  noise  tolerance  is  expected  for  this  IMMS  method.  At  ten  points  chosen  randomly  over  the  entire  image, 
good  intensity  variation  curves  are  recorded  by  driving  the  PZT  mirror  in  an  in-plane  DSPI  setup  where  the  setup  and 
specimen  are  all  on  a  vibration  insulated  table.  Numerical  simulations  and  some  real  noise  signals  are  used  to  test  the  noise 
tolerance  property  of  this  IMMS  method  at  these  ten  points. 


1.  Numerical  simulation 

A  numerical  Gaussian  white  noise  generator  whose  standard  deviation  can  be  controlled  is  used  to  generate  some  white  noise 
signals.  These  noise  signals  are  added  to  the  recorded  intensity  signals.  The  IMMS  method  is  employed  to  both  the  originally 
recorded  intensity  signals  and  the  noise-added  intensity  signals  to  determine  phase  values.  Results  are  compared  as  follows. 
Figure  1  is  a  recorded  good  intensity  variation  curve.  Figure  2  is  an  example  of  Gaussian  white  noise  whose  standard 
deviation  is  2.0.  The  RMS  signal  to  noise  ratio  is  6.3. 


Original  Wave 
10.6  z 


5.0- 

0.0 

-5.0- 


-9.9- 

0.0  20.0  40.0  60.0  80.0  100.0 


Figure  1 .  A  recorded  good  intensity  changing 
nur\/p 


Figure  2.  Gaussian  white  noise  signals. 


Figure  3  presents  the  fitted  curve  of  the  originally  recorded  intensity  wave.  Figure  4  shows  the  resulting  curve  from  IMMS  for 
the  noise-added  (Figure  2)  intensity  signals.  The  results  match  well. 
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Figure  3.  The  fitting  curve  of  the  original 
recorded  signals. 


Figure  4.  The  fitting  curve  of  noise  added 
intensity  signals. 


The  deviation  of  the  noise  signals  was  increased  gradually.  Figure  5  shows  the  fitting  result  for  the  same  recorded  intensity 
signals  as  in  Figure  1  after  a  noise  signal  whose  deviation  is  6.5  was  added.  Figure  5  shows  that  the  IMMS  method  cannot 
efficiently  process  such  a  noisy  signal.  Testing  through  a  range  of  deviation  of  the  white  noise  signals  shows  that  the  IMMS 
method  can  efficiently  remove  noise  having  only  small  deviation. 


The  same  noise  having  a  6.5  deviation  was  then  filtered  by  a  high  pass  filter,  which  means  only  high  frequency  noise  is  left  to 
be  added  to  the  original  signal.  Figure  6  shows  that  a  good  result  is  obtained.  This  indicates  that  high  frequency  noise  even 
with  large  amplitude  can  be  removed  by  the  IMMS  efficiently.  Repetition  of  the  simulation  tests  showed  that  the  IMMS  method 
has  very  good  immunity  to  both  high  frequency  noise  and  small  amplitude  low  frequency  noise.  Large  amplitude  low  frequency 
noise  can  damage  the  efficiency  of  the  noise  immunity  of  the  IMMS  method. 
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Figure  5.  The  fitting  curve  of  an  intensity  wave 
with  a  6.5  deviation  noise  added. 
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Figure  6.  The  fitting  curve  of  the  intensity 
wave  with  a  6.5  deviation  high  frequency 
noise  added. 


2.  The  effect  of  real  noise  signals 

At  various  points  on  the  specimen,  the  intensity  changes  induced  by  environmental  noise  (such  as  air  turbulence  and 
vibration)  ware  obtained  under  three  different  situations  by  the  DSPI  setup.  The  IMMS  method  was  employed  to  process  these 
noise-added  signals. 

A.  The  specimen  is  placed  on  a  vibration-insulated  table: 

A  noise  signal  example  is  shown  in  Figure  7.  This  noise  is  added  to  a  good  intensity  variation  curve  (Figure  8).  Figure  9  gives 
the  fitting  result  of  the  IMMS  method.  The  phase  value  obtained  from  the  IMMS  is  -3.14,  which  exactly  matches  the  phase 


obtained  from  the  original  signals.  In  this  example,  the  RMS  signal-to-noise  ratio  is  7.22  Rad.  On  this  table,  the  vibration  is 
minute.  So,  the  RMS-SNR  is  bigger  than  for  other  cases.  The  IMMS  works  very  well  in  this  situation,  as  expected. 
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Figure  7.  Noise  signals  under  A  condition. 


Figure  8.  Intensity  change  curve  with  noise 
added. 
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Figure  10.  Noise  signal. 


B.  The  specimen  is  put  on  the  same  table  without  air  support  vibration  isolation: 

Figure  1 0  is  an  example  of  noise  signals  under  this  situation.  The  original  intensity  curves  is  distorted  more  than  the  previous 
one  by  the  noise  (Figure  11).  Figure  1 2  shows  the  IMMS  result.  The  RMS-SNR  is  1 .81 .  The  phase  result  still  matches  exactly 
at  -2.29  Rad.  In  this  situation,  the  specimen  vibration  is  much  more  serious.  Therefore,  the  RMS-SNR  is  smaller  than  the 
previous  case,  around  1 .8.  Flowever,  the  IMMS  still  works  well. 
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Figure  1 1 .  Intensity  curve  with  noise  added. 


Figure  12.  The  IMMS  result. 


C.  Noise  recorded  with  the  specimen  placed  on  a  cart. 

In  this  case,  the  specimen  was  placed  on  an  ordinary  laboratory  cart  that  is  very  susceptible  to  environment  noise,  especially 
the  vibrations.  Noise  signals  were  recorded  at  various  points  and  much  bigger  variances  were  observed  for  these  noise 
signals,  as  expected  (see  Figure  13,  which  has  a  30.62  variance). 


Figure  13.  Noise  sample  obtained  under  case  C. 


In  this  situation,  an  alternative  calibration  was  used  to  measure  in-plane  displacement  of  a  plate.  The  recorded  noise  signals 
were  added  to  each  pixel  over  the  whole  image  during  displacement  measurement  by  the  IMMS  algorithm.  The  same 
displacement  calibration  set  up  as  in  reference  [10]  was  used  to  introduce  known  in-plane  displacement.  Table  1  gives  some 
examples  of  these  test  results.  It  indicates  that  the  IMMS  worked  well  even  under  these  high  noise  conditions. 


Table  1 .  Some  examples  of  displacement  measurements 


Variance  of  noise 

Induced  known 
displacement  (pm) 

Measured  displacement  by 
the  IMMS  (pm) 

Error  (%) 

30.62 

1.864 

1.85 

-0.75107 

37.75 

1.682 

1.75 

4.042806 

30.23 

2.704 

2.82 

4.289941 

All  the  above  simulations  were  repeated  tens  of  times  at  different  dates  and  times.  The  IMMS  method  proved  to  have  a  good 
ability  to  reduce  certain  specific  noise  effects.  Especially,  the  IMMS  method  reliably  removes  high  frequency  noise  and  small 
amplitude  low  frequency  noise.  The  noise  immunity  ability  of  the  newly  developed  IMMS  method  makes  the  DSPI  technique 
more  practical. 


4.  Application  of  IMMS  to  PZT  calibration 

By  applying  the  phase  determination  computation  presented  in  section  2  at  all  points  along  the  entire  intensity  curve,  the  plot 
of  phase  values  versus  driving  voltages  is  obtained  with  2n  ambiguity  as  shown  in  Figure  14.  The  phase  curve  is  obtained  by 
performing  phase  unwrapping  with  respect  to  the  zero  driving  voltage  point.  The  phase  steps  induced  by  the  driving  voltages 
are  determined  (Figure  15).  From  Figure  15,  the  driving  voltage  required  by  the  PZT  to  get  an  expected  phase  step  at  this 
point  can  be  obtained  by  interpolation  of  the  data.  As  well,  the  non-linearity  of  the  phase  shifting  at  this  point  is  displayed  on 
the  curve.  To  calibrate  a  specific  phase  step,  i.e.  n/2  phase  step,  the  above  calibration  can  performed  at  a  number  of  pixel 
points  over  the  entire  image.  An  average  driving  voltage  to  get  that  phase  step  can  be  obtained. 


Figure  14.  Phase  Step  with  2k  ambiguity. 


Figure  15.  Phase  Step  VS  Driving  Voltage. 


Tilt  and  non-uniformity  of  the  phase  shifter  movement  are  two  other  important  concerns  that  significantly  affect  the  accuracy  of 
phase  shifting.  In  this  study,  at  a  given  driving  voltage  to  the  PZT,  phase  shifting  data  at  points  along  two  perpendicular  edges 
of  the  image  were  computed  and  plotted.  The  calculation  procedure  described  above  was  repeated  at  each  point  on  these  two 
edges  to  get  the  phase  shifting  curves.  Figure  16  shows  the  tilt  and  non-uniformity  of  phase  shift  along  the  x-axis  edge,  and 
Figure  17  is  for  the  y-axis  edge.  The  phase  step  values  corresponding  to  the  calibrated  n/2  driving  voltage  were  combined  into 
a  3-D  graph,  from  which  the  tilt  and  non-uniformity  of  the  phase  shifter  can  be  assessed  visually  and  directly.  Figure  1 8  is  an 
example  of  such  a  3-D  graph. 
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Figure  16.  Tilt  and  uniformity  in  x  direction. 
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Figure  1 7.  Tilt  and  uniformity  in  y  direction. 


Figure  18.  A  3-D  plot  of  the  obtained  phase  shifting  results. 
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The  method  described  above  has  been  used  in  our  lab  to  calibrate  and  inspect  the  phase  shifter  used  in  a  DSPI  setup.  Figure 
19  represents  a  recorded  intensity  wave.  After  filtering  and  curve  fitting,  a  smooth  waveform  was  obtained  as  shown  in  Figure 
20. 


Figure  19.  Recorded  intensity  waveform 


Figure  20.  Fitting  curve  obtained  by  IMMS 


The  calibration  result  is  as  shown  in  Figure  15.  For  this  specific  experimental  setup,  the  driving  voltage  to  get  a  n/2  phase  step 
is  2.09V.  Figure  1 5  shows  that  the  linearity  of  this  phase  shifter  is  good  as  well.  Figures  16,17  and  1 8  indicate  that  the  tilt 
and  non-uniformity  of  this  PZT  phase  shifter  are  very  small. 

From  the  results  of  hundreds  of  tests,  the  average  driving  voltage  to  get  a  n/2  phase  step  is  determined  at  2.13V  for  this 
apparatus.  The  standard  deviation  is  0.104,  which  implies  very  good  stability  for  this  setup. 

5.  Conclusions  and  Comments 

In  this  paper,  the  noise  tolerance  of  IMMS  method  is  discussed.  It  is  shown  that  the  IMMS  method  has  a  good  tolerance  to 
certain  noise  signals,  including  high-frequency  noise  and  small-amplitude  low-frequency  noise.  Thus,  the  IMMS  technique 
facilitates  application  of  interferometric  techniques,  including  DSPI/ESPI,  in  factory-type  environments. 

A  phase  shifting  calibration  algorithm  has  been  derived  from  the  IMMS  method  also,  and  its  applications  to  a  PZT  phase 
shifter  were  successfully  demonstrated  in  some  experiments.  This  algorithm  shows  several  practical  advantages.  First,  it 
eliminates  the  drawbacks  caused  by  the  common  assumptions  of  linearity  and  uniformity  of  the  phase  shifter,  thus  efficiently 
minimizing  the  movement  error  of  the  PZT.  Second,  it  can  be  used  to  inspect  the  working  reliability  of  the  phase  shifter  by 
showing  simultaneously  the  non-linearity,  non-uniformity,  and  tilt  of  the  phase  shifter. 

Overall,  this  new  method  is  a  convenient  and  ideal  way  to  calibrate  the  phase  shifter,  as  well  as  to  inspect  the  quality  of  the 
phase  shifter  and  the  reliability  of  the  experimental  setup.  Whenever  the  experimental  setup  is  adjusted,  such  as  by  changing 
the  illumination  angle  or  the  path  length  of  beam,  this  algorithm  can  implement  calibration  and  inspection  online. 
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